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Abstract. We present new developments on the Cosmic-Ray driven, galactic dynamo, mod- 
eled by means of direct, resistive CR-MHD simulations, performed with ZEUS and PIERNIK 
codes. The dynamo action, leading to the amplification of large-scale galactic magnetic fields 
on galactic rotation timescales, appears as a result of galactic differential rotation, buoyancy 
of the cosmic ray component and resistive dissipation of small-scale turbulent magnetic fields. 
Our new results include demonstration of the global-galactic dynamo action driven by Cosmic 
Rays supplied in supernova remnants. An essential outcome of the new series of global galactic 
dynamo models is the equipartition of the gas turbulent energy with magnetic field energy and 
cosmic ray energy, in saturated states of the dynamo on large galactic scales. 
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1. Introduction 

Numerous processes have been proposed to explain initial magnetic fields in early 
galaxies. However all these processes, like phase transitions in the early universe or 
Biermann battery in protogalactic objects provide only very week magnetic fields of 
the order of B ~ 10~ 20 G at the beginning of galactic evolution. On the other hand, 
Rees (j!987j) proposed that initial magnetic fields in galaxies might have been generated 
in first stars and then scattered in the interstellar medium (ISM) by SN explosions, 
and subsequently amplified in plerionic (Crab-type) supernova remnants (SNRs). The 
mean magnetic field on the galactic scale has been estimated to be of the order of 10 -9 G. 
Radio observations (see Beck Be ck 2009} this volume) indicate that typical contemporary 
magnetic fields in spiral galaxies are of the order of few up to few tens of /zG, which means 
that magnetic fields have undergone amplification by at least four orders of magnitudes 
within the galactic lifetime. Therefore, a model for efficient magnetic field amplification 
during galactic evolution is necessary. The standard model of magnetic field amplification 
in disk galaxies is based on the theory of turbulent mean field dynamo (see Widrow 2002 
for a recent review). 

The mean-field dynamo theory has been successful in explaining magnetic fields in var- 
ious astrophysical objects, however, the classical, kinematic dynamo seems to be rather 
slow in galaxies. Magnetic field amplification timescale ^dynamo ~ (0.5 1) x 10 9 yrs is 
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too long to explain ~ 1/j,G magnetic fields in galaxies which are only few 10 9 yrs old 
(z ~ 1 3) (Wolfe, Lanzetta & Oren [19921 Kronberg et al. 12008)1 . These circumstances 
seem to indicate a need of an alternative approach. We suggest that direct numerical 
MHD modeling of ISM dynamics, through considerations of detailed physical processes, 
involving major ISM components: gas, magnetic fields and cosmic rays, is necessary to 
follow magnetic history of galaxies. 

2. Cosmic Rays in the interstellar medium 

The dynamical role of CRs was recognized for the first time by Parker (1966), who 
noticed that a vertically stratified ISM, consisting of thermal gas magnetic field and CRs 
is unstable due to buoyancy of the weightless components: magnetic fields and CRs. 
The CR component appears to be an important ingredient of this process. According 
to the diffusive shock acceleration models CRs are continuously supplied to ISM by SN 
remnants. Therefore, strong buoyancy effects due to CRs are unavoidable. 

Theories of diffusive shock acceleration predict that about 10 % of the ~ 10 51 erg of 
the SN II explosion energy is converted to the CR energy (see e.g. Jones et al 1998.1 and 
references therein). Observational data indicate that gas, magnetic fields and CRs appear 
in approximate energetic equipartition, which means that all the three components are 
dynamically coupled. Moreover, numerical experiments by Giaccalone & Jokipii (1999) 
indicate that CRs diffuse anisotropically in the ISM, along the mean magnetic field 
direction. 

3. Cosmic Ray transport 

To incorporate CR propagation in MHD considerations we use the diffusion - advection 
equation (e.g. Schlickeiser & Lerche 1985) 



where the diffusion term in written in the tensorial form to account for anisotropic 
diffusivity of CRs. The sources of CRs on the rhs. of Eqn. 13.11 correspond to the CR 
production in Supernova remnants. 

The dynamics of ISM including the energetics of CRs can be now described by the CR 
transport Eqn. Q3.ip and the system of MHD equations. We note that in the presence 
of cosmic rays additional source term: — VPcr should be included in the gas equation 
of motion (see e.g. Berezinski et al. 1990), in order to incorporate the effects of CRs on 
gas dynamics. The cosmic ray diffusion-ad vection Eqn. (|3.ip has been supplemented to 
the MHD algorithm of Zeus-3D MHD code by Hanasz & Lesch (j2003ap with the aim of 
studying ISM physical processes, e.g. Parker instability, in which CRs play a dynamically 
important role. 

4. CR-driven dynamo 

Referring back to the issue of galactic magnetic fields, we are primarily interested in 
processes leading to the amplification of large scale magnetic fields. The original idea of 




cr 



+ V(e cr u) = -p cr V • v + V(KVe cr ) 
+ CR sources (SN remnants) 



(3.1) 
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CR-driven dynamo has been raised by Parker (1992). Our first local (shearing-box) CR- 
driven dynamo numerical experiments (Hanasz et al. 2004^ 2006; Otmianowska-Mazur et 
al. 2007) rely on the following ingredients: (1.) the cosmic ray component described by the 
diffusion-advection transport equation (|3.ip . including localized sources of cosmic rays 
— supernovae remnants, exploding randomly in the disk volume, while SN shocks and 
thermal effects are neglected (see Gressel et al. 2008a J 2008b) for a recent dynamo model 
relying on the thermal energy output from SNe); (2.) resistivity of the ISM (see Hanasz, 
Otmianowska-Mazur and Lesch 2002), leading to topological evolution of magnetic fields; 
(3) shearing boundary conditions (Hawley, Gammie and Balbus 1995) together with 
Coriolis and tidal forces, aimed at modeling of differentially rotating disks in the local 
approximation; (4) realistic vertical disk gravity and rotation following the model of ISM 
in the Milky Way by Ferriere (j!998j) . 

In the most recent paper (Hanasz et al. 2009) we present a parameter study of the 
CR-driven dynamo model by examining the dependence of magnetic field amplification 
on magnetic diffusivity, supernovae rate determining the CR injection rate, temporal 
modulation of SN activity, grid resolution, and CR diffusion coefficients. We find the 
dominating influence of magnetic diffusivity (treated as a free parameter), among other 
parameters, (Fig. [1]) on the efficiency of magnetic field amplification. 
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Figure 1. Time evolution of azimuthal magnetic flux and total magnetic energy for different 
values of magnetic diffusivity in simulation series A. The curves represent respectively cases of 
77 = (Al), 77 = 1 (A2), 77 = 10 (A3), 77 = 100 (A4) and 77 = 1000 (A5) in units pc 2 Myr -1 . 

The fastest magnetic field amplification, observed in the experiments, coincides with 
the magnetic diffusivity comparable to n ~ \r]turb{obs), where r] tU rb{obs) ~ f 100 pc x 
10 km s _1 ~ 10 26 cm s -1 . The question of which physical process can be responsible 
for that large magnetic diffusivity coefficients points out towards the current studies on 
magnetic reconnection (Lazarian et al. I2004[ Cassak et al. 2006) 

To investigate observational properties of current dynamo models we perform a series 
of numerical simulations, placing the local shearing-boxes at different galactocentric 
radii. We replicate the contents of local cubes into rings, and then combine rings to 
build up a disk (Otmianowska-Mazur et al. 2009). A synthetic polarization radio-map 
of synchrotron emission constructed, on the base of simulated magnetic fields and CR 
distribution is shown in Fig. [U 



122 



Hanasz et al. 




14 54 06 04 02 00 S3 SS 56 54 52 50 J4 54 06 04 02 00 53 5S 56 54 52 50 
RIGHT ASCENSION (J2000) L{|( HIT .\S< I NSh >N .12 



Figure 2. Synthetic radio-map of polarized radio emission for our model at the inclination 
i = 80° and a map of real galaxy NGC 5775 

A polarized emission map together with superposed polarization vectors of synchrotron 
emission of an edge-on galaxy NGC 5775 is shown for comparison. As it is apparent 
the X-shaped structures of magnetic field, present in the real galaxy, form also in the 
synthetic radio image. The presented results demonstrate that the CR-MHD simulations 
provide a way of verification of CR-dynamo models against the observational data. 

To summarize the overall outcome of the local simulations of CR-driven dynamo, we 
note that the model provides efficient amplification of large-scale magnetic fields, which 
attain, in the saturated state, equipartition with turbulent gas motions. However, we note 
also a drawback of the local models, since the CR energy density remains more than an 
order of magnitude larger than gas and magnetic energy densities. We recognize that 
improvements of the current models are necessary, since the possible reason of the excess 
of cosmic rays may result from the effect of trapping of cosmic rays by predominantly 
horizontal magnetic field in the computational box which is periodic in horizontal direc- 
tions. Therefore, global galactic disk simulations should be performed to enable escape 
of CRs along the predominantly horizontal magnetic field in the galactic plane, in order 
to avoid the CR excess. 

5. Global disk simulations 

Recently we started a new series of global simulations of CR-driven dynamo with the 
aid of PIERNIK MHD code (see Hanasz et al. I2009a[ I2009b[ I2009c[ I2009dp . which is a 
grid-based MPI parallelized, resistive MHD code based on the Relaxing TVD (RTVD) 
scheme by Jin & Xin (1995) and Pen et al. (|2003p . The original scheme is extended to 
deal with the diffusive CR component (see Hanasz & Lesch 2003a). 

In Fig. [3] we show results of one of the first semi-global simulations of CR-driven 
dynamo. The simulations have been performed for a quarter of galactic disk resembling 
Milky Way, with parallel CR diffusion coefficient K\\ = 3 x 10 28 cm 2 s _1 , in a computa- 
tional domain 20 kpc x 20 kpc x 20 kpc, and resolution of 400 x 400 x 160 grid cells, in 
16 MPI blocks. 

The semi-global simulations confirm the fast amplification rate of the mean magnetic 
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Figure 3. Horizontal and vertical slices through the computational box for a simulation of 
CR-driven dynamo in a semi-global galactic disk. The color maps reflect: 1: gas density, 2: CR 
energy density, 3: magnetic field strength, 4. the ratio of magnetic to CR energy densities and 5. 
synchrotron emissivity for a saturated state of the dynamo. The gray-scales in all panels, except 

4. are chosen arbitrarily. The darker shades correspond to larger values of displayed quantities. 
The greyscale of panel 4. is such that black areas represent values of e mag /e cr ~ 1. 

field: growth time of the mean magnetic flux is approximately 210 Myr in the whole 
galaxy, which is close to the galactic rotation period near the orbit of Sun. As expected, 
the global simulations introduce a significant improvement with respect to the former 
shearing-box simulations. It is apparent in Fig. [3] that the ratio of magnetic to cosmic 
ray energy is close to one in the synchrotron emitting volume (compare panels 4. and 

5. ), while over-equipartition of CR energy with respect to magnetic energy still holds 
outside the disk. 

The latest development of our CR-dynamo model is a fully global galactic disk simula- 
tion (see the complementary paper by Hanasz, et al. this volume), where we demonstrate 
that dipolar magnetic fields supplied on small SN-remnant scales, can be amplified expo- 
nentially by the CR-driven dynamo to the present equipartition values, and transformed 
simultaneously to large galactic-scales. 



6. Conclusions 

We have shown that the CR contribution to the dynamics of ISM, studied by means of 
CR-MHD simulations, in both local and global scales, leads to a very efficient magnetic 
field amplification in galactic disks. 

The cosmic ray driven dynamo amplifies efficiently galactic magnetic fields on a ti- 
mescale of galactic rotation. The resulting growth of the large-scale magnetic field by 4 
orders of magnitude within 2 Gyr (Hanasz at al. 12004)) . is fast enough to expect ~ 1/zG 
magnetic field in galaxies at z ~ 1 3. 

We point out the advantage of global disk models, which on the contrary to shearing- 
box models, provide CR energy equipartition with magnetic field in the synchrotron 
emitting part of the disks. 
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